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ABSTRACT 

Combined Spitzer, Chandra, XMM-Newton, and VIA observations of the giant radio galaxy NGC 1316 
(Fornax A) show a radio jet and X-ray cavities from AGN outbursts most likely triggered by a merger with 
a late-type galaxy at least 0.4 Gyr ago. We detect a weak nucleus with an SED typical of a low-luminosity 
AGN with a bolometric luminosity of 2.4 x 10 42 ergs s" 1 . We examine the Spitzer IRAC and MIPS images of 
NGC 1316. We find that the dust emission is strongest in regions with little or no radio emission and that the 
particularly large infrared luminosity relative to the galaxy's K-band luminosity implies an external origin for 
the dust. The inferred dust mass implies that the merger spiral galaxy had a stellar mass of 1 - 6 x 10 10 M Q and 
a gas mass of 2-4 x 10 9 M Q . X-ray cavities in the Chandra and XMM-Newton images likely result from the 
expansion of relativistic plasma ejected by the AGN. The soft (0.5-2.0 keV) Chandra images show a small 
~ 15" (1.6 kpc) cavity coincident with the radio jet, while the XMM-Newton image shows two large X-ray 
cavities lying 320" (34.8 kpc) east and west of the nucleus, each approximately 230" (25 kpc) in radius. Current 
radio observations do not show emission within these cavities. The radio lobes lie at radii of 14.'3 (93.3 kpc) 
and 15 '6 (101 kpc), more distant from the nucleus than the detected X-ray cavities. The relative morphology of 
the large scale 1 .4 GHz and X-ray emission suggests they were products of two distinct outbursts, an earlier one 
creating the radio lobes and a later one producing the X-ray cavities. Alternatively, if a single outburst created 
both the X-ray cavities and the radio lobes, this would require that the radio morphology is not fully defined by 
the 1 .4 GHz emission. For the more likely two outburst scenario, we use the buoyancy rise times to estimate an 
age for the more recent outburst that created them of 0. 1 Gyr and the PV work done by the expanding plasma to 
create the X-ray cavities to estimate the outburst's energy to be 10 58 ergs. The present size and location of the 
radio lobes implies that the outburst that created them happened ~ 0.4 Gyr ago and released ~5 x 10 58 ergs. 

Subject headings: galaxies: active — galaxies: individual (NGC 1316) — galaxies: structure — infrared: 
galaxies — radio continuum: galaxies — X-rays: galaxies 



1. INTRODUCTION 

NGC 1316 (Fornax A) is one of the near est and brightest ra - 
dio galaxies with radio lobes spanning 33' (lEkers et al.ll983h . 
It lies on the outskirts of the Fornax cluster at a distance of 
22.7±1.8 Mpc, based on a distance modulus of 31.66±0.17 
dTonrv et aljfeoOll) Early optical observations of NGC 1 3 1 6 
( Evansi 1949b r e veale d dust lanes in the nuclear region, lead- 
ing IShklovskiil (119631) to hypothesize that radio lobes might 
be powered by accretion of interstellar gas onto the nucle us. 
More extensive optical observations led Schweizer ( 1980) to 
classify NGC 1316 as a D-type galaxy dMorganlll958l) with 
an elliptical-like spheroid embedded in a large envelope. 

Schweizer further suggested that NGC 1316's disturbed 
morphology may be due to one or mor e low-mass gas-rich 
merge rs, occurring over the last 2 Gyr. iMackie & Fabbiano 
( 1998) concluded, from more recent optical observations, in 
combination with ROSAT data, that NGC 1316 had undergone 
either a major merger more than 1 Gyr ago or a merger with 
a low-mass gas-rich galaxy about ~ 0.5 Gyr ago. NGC 1316 
exhibits other signs of at least one merger, including loops of 
Ha filaments resembling tidal tails with projected lengths as 
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large as 10' (65.4 kpc) and a gas disk rotating much faster 
than the stellar spher oid and at an angle to it, indicating a 
likely external origin dXilouris et al.ll2004l ISchweizerl 1980). 
Goud frooii et alj (1200 lb used globular clusters to date any ma- 
jor mergers to between 1 .5 and 4 Gyr ago, but the age and type 
of merger (or mergers) still remain uncertain. 

Multiwavelength observations can set constraints on the 
merger event. Radio emission, from the radio jet and lobes, 
is expected to be powered by accretion onto the central super- 
massive black hole (SMBH), which can be enhanced by the 
infall of material from a merger. The expanding radio plasma 
can create cavities in the surrounding hot gas, which are seen 
as decrements in the X-ray emission. By m easuring the PV 
work done by the expand ing radio plasma (McNamar a et al.l 
120001 IChurazov et al. 2002), we can constrain the energy pro- 
duced by the SMBH. Also, since the X-ray cavities rise buoy- 
antly, their distance from the nucleus constrains the age of the 
outburst. We can set mass-related constraints on the merger 
galaxy by examining NGC 1316 for dust. As an early-type 
galaxy, NGC 1316 is expected to be dust poor. Mid-infrared 
observations permit us to measure the amount of dust present, 
which, if higher than expected for a galaxy of its size and type, 
indicates an external origin for the dust. 

In this paper, we report our analysis of the surface bright- 
ness distribution of NGC 1316 in the mid-IR with the Infrared 
Array Camera (IRAC) and the Multiband Imaging Photome- 
ter for Spitzer (MIPS) of the Spitzer Space Telescope and the 
resulting determination of the warm dust morphology. We 
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describe our analysis of the X-ray emission imaged both by 
Chandra and XMM-Newton. Observations and data reduc- 
tion are described in §2, and imaging and modeling results 
are given in §3. We compare the features seen in the infrared 
and X-ray with the radio emission in the nuclear, inner jet, and 
extended emission regions in §4 and examine the constraints 
these observations place on the mass of the merger progen- 
itor, the outburst and merger ages, and the outburst energies 
in §5. Finally, we summarize the results in §6. Images have 
north to the top and east to the left. Angles are given counter- 
clockwise from west, unless otherwise stated. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Spitzer Observations 

IRAC dFazio et all 12004b observations of NGC 1316 were 
obtained on 2004 July 19 and 22 as part of the SIRTF Nearby 
Galaxy Survey (SINGS) Legacy program (Kennicutt et al. 
2003; Spitzer PID 159) in all four bands. The two visits to 
NGC 1316 consisted of similar 4 x 30 s integrations covering 
the galaxy, its companion NGC 1317, and the nearby field. 

For our analysis, we retrieved the Basic Calibrated 
Data (BCD) version S14.0 pipeline products from the 
Spitzer archive. We corrected the IRAC BCD frames for 
residual images arising from prior, unrelated observations 
of bright sources by making object-masked, median-stacked 
coadds of all science frames not containing significant ex- 
tended emission from NGC 1316orNGC 1317. These single- 
visit coadds were then subtracted from the individual BCDs of 
the corresponding visit to remove the residual images. For the 
8.0 /im observation, light scattered from the galaxy nucleus 
along the detector array rows was fit and subtracted from th e 
BCD frames using custom software as in Ash bv et al.l d2009). 
We then coadded the modified 30 s BCD frame s from both ob- 
serva tions using version 4.1.2 of IRACproc (Schust er et al.l 
2006) to mosaics having l."2 pixels, i.e., the native IRAC 
pixel size. 

The SIN GS program also o btained MIPS observations of 
NGC 1316 (iRieke etal .1 12004) on 2004 December 5 and 7. 
We obtained a combined 24 pm mosaic and coverage map of 
NGC 1316 from the SINGS websitsfl 

A 12' x 12' region centered on the galaxy was selected 
for analysis in each IRAC and MIPS mosaic. The resulting 
3.6 /im image, which is very similar to the 4.5 pm image, is 
presented in FigureQJ. The corresponding images for 5.8 pm, 
8.0 /im, and 24 /im are shown in the top row of Figure [2] Un- 
certainty maps were created by adding pixel-pixel rms noise 
and shot noise in quadrature. 

2.2. Chandra Observations 

NGC 1 3 1 6 was observed for 30 ks on 200 1 April 1 7 (ObsID 
2022) with the Chandra Advanced CCD Imaging Spectrom - 
eter Spectroscopy Array (ACIS-S; Weiss kopf et al.l [2 000). 
This o bservation was previously studied by iKim & Fabbianol 
(2003). We reprocessed this observation applying the latest 
CTI and time-dependent gain calibrations (see Vikhlinin et 
al. 2005 for more details). We applied the standard filter- 
ing by grade, excluded bad/hot pixels and columns, and re- 
moved cosmic ray 'afterglows'. We removed time intervals 
with background flaring, resulting in an effective exposure 
time of 23.853 ks. The background files (see Markevitch et 
al. 2000 for details) were processed in exactly the same man- 
ner as the observations. 

4 http://data.spitzer.caltech.edu/popular/sings 



2.3. XMM-Newton Observations 

XMM-Newton observed NGC 1316 for 107.3 ks on 2005 
August 11-12 (ObsID 0302780101). He re we examine the 
Metal Oxide Semi-conductor CCD (MOS: Uansen et al.l200lh 
observation. We filtered the data to remove periods of back- 
ground flaring, resulting in a reduced exposure time of 62.0 ks 
for MOS 1 and 56.3 ks for MOS2. The events were further fil- 
tered to retain only events with energies between 0.5 and 7.0 
keV and patterns between and 12. 

2.4. Radio, Optical, and CO Observations 

NGC 1316 was observed with the Very Large Array (VIA) 
at 4.89 GHz in the AB array configuration on 2002 June 1. 
We used the Astronomical Image Processing System (AIPS; 
version 31DEC09) package to generate a map with a res- 
olution of 1"38 x 1"02 (Figure [3};). We also obtained a 
20 cm VIA map from the NAS A/IPAC Extragalactic Database 
(NED) dFomalont et al.|[l989h . 

The Hubble Space Telescope Advanced Camera for Sur- 
veys (ACS) observed NGC 1316 on 2003 March 04 through 
the F555W filter for 6.98 ks (PropID 9409). The WFPC2 in- 
strument on Hubble observed it through the F814W filter on 
1996 April 07 for 1.86 ks (PropID 5990). I2 CO(2-l) inten- 
sities at 230 GHz were obtained from C. Horellou and J. H. 
Black, who observed NGC 1316 with the 15 m Swedish-ESO 
Submillimeter Telescope (SEST) in 1999 and 2001 with a res- 
olution of 22" dHorellouetal.1120011) . 

3. DATA ANALYSIS 

3.1. Spitzer Analysis 

To measure dust and nuclear emission, we first had to re- 
move the stellar component from the Spitzer images. To ac- 
complish this, we modeled the 3.6 pm emission, where the 
flux from the stars is the greatest of all four IR AC bands, with 
a two-dimensional Sersic model (Sersic 119681) for the stellar 
contribution. To determine the galaxy model, as well as emis- 
sion from the nucle us, we used the software package GALFIT 
dPeng et al. 1120021) . which is a parametric surface brightness 
fitting code using \ 2 minimization. We iteratively determined 
the center for the Sersic profiles based on the 3.6 pm emis- 
sion and the position of the central point source based on the 
8 pm emission. After confirming that the 4.5 pm emission 
results in similar parameters, as expected, since this band is 
also dominated by stellar emission, we held the Sersic index, 
effective radius, axis ratio, position angle, and central posi- 
tion fixed at the 3.6 /im fitted values for subsequent fits to the 
4.5 pm, 5.8 pm, 8.0 /im, and 24 pm images. The remaining 
free parameters were the normalizations of the Sersic model, 
central point source, and the foreground star 35" southeast of 
the nucleus, as well as a sky model, consisting of a constant 
offset and gradients in the two array directions. The point re- 
sponse function (PRF) in each band was input to GALFIT in 
order to fit the emission from the nucleus and foreground star. 
Masks were used to exclude foreground stars and the bulk of 
the emission from the neighboring galaxy NGC 1317. The 
fitted model parameters are given in Table Q] The fitted Ser- 
sic index of 6.07 ±0.10 agrees well with indices of 5.8 and 5.9 
determined by Cote et al. (2007) from Hubble ACS images of 
NGC 1316. 

Figure QJ shows the 3.6 fim mosaic. Figure QJ shows the 
best fitting model obtained from GALFIT of a central point 
source, foreground point source, and a Sersic model. Since 
the residual emission, shown in Figure QJ, is on the order of 
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FIG. 1 . — The input image used in GALFIT when fitting 3.6 fim (a), the final model combining a Sersic profile, a point source, and sky (b), and the fit residuals 
(c). The 4.5 fim images are similar to the 3.6 fim images, and are therefore not shown here. Each image is 3'0 X 3'0 with north to the top and east to the left, as 
are all images throughout this paper unless stated otherwise. The color scale for all the panels is inverted so bright regions in the residual image (c) are regions 
where the model is over-subtracted. Panels a and b have logarithmic color scales, while panel c has a linear scale. The dark point source 35" southeast of the 
galaxy's nucleus is a foreground star, which we also fit. The scale bar corresponds to 30" (3.3 kpc). 
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FIG. 2. — The top row shows the input image used in GALFIT when fitting the 5.8 fim (left), 8.0 fim (middle), and 24 fim (right) images. The bottom row are 
those same images, once a Sersic model and the foreground star have been subtracted. The image sizes are the same as in Figure[T] The color scale is inverted, 
and the top row has a logarithmic color scale, while the bottom row has a lineal' color scale. The dark regions in the residual images are excess emission, most 
likely due to dust. The scale bar corresponds to 30" (3.3 kpc). 
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the noise (~ 3%), outside the nuclear region, our Sersic model 
is a good stellar emission model for this galaxy. Since the 
4.5 /im image and residuals are similar to those at 3.6 /im, 
they are not shown. 

Figure [2] presents the 5.8 /im, 8.0 /im, and 24 /im im- 
ages (top row) and residuals after the modeled Sersic profile 
has been subtracted (bottom row). The 8.0 /im image and 
contours of its n onstellar emission were previously shown in 
Temi et alJ d2005h . The non-axially symmetric, non-stellar 
component is visible in the top row (a-c) and quite striking 
in the bottom row (d-f). As we discuss in §4 and §5.1, most 
of this component is due to dust. In addition to the nucleus, 
there is a region approximately 11 "8 (1.3 kpc) in radius of 
dust emission 29" (3.1 kpc) southeast of the nucleus. There is 
also dust emission extending 44" (4.8 kpc) towards the north- 
west, ending in a clumpy arc that extends over a —90° angle 
and containing two knots 21 "8 (2.4 kpc) and 32."0 (3.5 kpc) 
from the nucleus. 

3.2. Chandra Analysis 

We analyzed Chandra observations of NGC 1316 to deter- 
mine the brightness and morphology of emission due to hot 
gas and a central nuclear source. We created an image of the 
soft X-ray emission between 0.5 and 2 keV using the ACIS- 
S3 CCD. We re moved point sources, whi ch we detected using 
WAVDETECT dFreeman & Vinavll2002l) . We then created an 
exposure map that accounts for all position dependent, but en- 
ergy independent, efficiency variations across the focal plane 
(e.g., overall chip geometry, dead pixels or rows, variation of 
telescope pointing direction). Finally, we made a flat-fielded 
image by subtracting both the blank-field background and the 
read-out background and then dividing by the exposure map. 
We also created an image of the hard X-ray emission be- 
tween 2 and 7 keV. The resulting images, shown in Figures 
01 and[3f, have 1" and 0."5 pixels respectively and have been 
smoothed with a 5" Gaussian. 

3.3. XMM-Newton Analysis 

We analyzed the XMM-Newton image to examine features 
of the extended X-ray emission which were not in the smaller 
Chandra field of view. We created new exposure maps and 
images, to remove the contribution from the MOS 1 CCD with 
higher background, and we removed point sources. The ex- 
tended X-ray emission, described in §4.3, includes a pair of 
cavities 320" (34.8 kpc) to the west and southeast of the nu- 
cleus. We reprojected blank sky background files for each 
of the MOS CCDs^ to the coordinates of our event file and 
scaled them appropriately, prior to subtracting the combined 
background file from the image. We then divided by the expo- 
sure map to create an exposure-corrected image with 4" pixels 
(Figure 0J. 

4. MULTIWAVELENGTH COMPARISON OF FEATURES 

Figure [3] shows NGC 1316 in visible (555 nm), non-stellar 
infrared (8.0 /im), radio (4.89 GHz), Chandra soft (0.5- 
2.0 keV), XMM-Newton soft (0.5-2.0 ke V), and Chandra hard 
(2.0-7.0 keV) emission. For the 8.0 /im image, we subtracted 
the stellar emission, as modeled by a Sersic profile. In Figure 
|5] we show the relative locations of the X-ray, IR, radio, and 
CO emission. Below, we discuss the nucleus, the jet, and the 
extended non-stellar emission. 

5 http://xmm.vilspa.esa.es/extemal/xmm_sw_cal/ 
background/blank_sky. shtml 



4.1. Nuclear Emission 

The nucleus is detected in the radio, UV, X-ray, and IR 
bands. We obtained radio and UV nuc l ear fluxes from the lit - 
erature (Geldz ahler & Fomalontll 19841 iFabbiano et al.lll994l) . 
We measured the X-ray spectrum in a 1 " radius circle around 
the nucleus, using a l"-2" annulus to subtract the background 
thermal emission. We fit the nuclear spectra with an absorbed 
power law with both a variable and a n held constant a t 
the Galactic value of 2.4 x 10 20 cm -2 dKalberla et al.ll2005h . 
The variable n# fit had a similar power law index to the 
constant n# fit and a column density not significantly larger 
than the Galactic absorption. Therefore, we used the constant 
iiH spectral fit whose best-fit power law index was 2.09_g 4 g 
(90% confidence), which is within the range typically ob- 
served for active galactic nuclei (AGN). We converted mea- 
sured soft (0.5-2.0 keV) and hard (2.0-7.0 keV) X-ray fluxes 
to luminosities of 6.5 ± 1.8 x 10 38 ergs s" 1 and 5.7±2.1 x 
10 38 ergs s . We measured the X-ray nuclear broadband 
(0.3-8 keV) flux to be 2.4 ±0.7 x 10" 14 ergs s" 1 cm" 2 . Our 
flux corresponds to a broadband (0.3-8.0 keV) luminosity of 

1.5 ±0.5 x 10 39 ergs s . We also measured the nuclear flux 
in the F814W Hubble image in a 0." 15 radius aperture to be 
7.8 ±0.2 x 10~ 13 ergs s -1 cm -2 , which we consider an upper 
limit because of the stellar contribution. Since the F555 W im- 
age is saturated in the central region, we do not measure the 
nuclear flux through this filter. 

The nucleus is detected in all five Spitzer bands. Fluxes and 
luminosities, vL v , in each photometric band from the GAL- 
FIT model are given in Table Q] We note that the nuclear 
fluxes were derived from fitting, but that the point source 
is not directly visible in the IRAC images. To confirm that 
the galaxy does contain a point-like nucleus, we used the 

3.6 //m emission as a stellar model, which was then color cor- 
rected and scaled to correct for differences in the apertures 
and zero point magnitudes between the IRAC bands. The 
resulting non-stellar images are very similar to those in the 
bottom row of Figure [2] The IR color of NGC 13 16's nu- 
cleus falls outside, but within 3a, of the region defined by 
the mid-IR AGN color selection criteria o flStern et alJ d2005l) 
and within the selection region defined by lLacy et alJ (12004). 
We note however that these IRAC selection criteria were de- 
veloped for Seyfert galaxies and quasars that are significantly 
more luminous than NGC 1316. Indeed, the spectral energy 
distribution (SED) of the NGC 1316 AGN (Figure© is sim- 
ilarto those of other low luminosity AGN (LLAGN) dHol 
1999) in that it lacks the big blue bump found in powerful, 
optically bright AGN and instead appears to show only a sin- 
gle big red bump and has a larger r adio to o ptical ratio than 
that of the higher luminosity AGN dHoll2008l) . The compar- 
ison s pectra plotted in this figure are from Fig ure 7 of iHo l 
(2008) where the purple squares are the LLAGN. Smi th et al] 
(120071) cited the AGN of NGC 1316 as having the typical 
peculiar PAH spectrum of LLAGN, which has low ratios of 
L(7.7/im)/L(11.3 /im). 

We used the supermassive black hole mass of 1.5 x 
10 8 M dNowak et alj|2008l) to determine the Eddington lu- 
minosity for the AGN to be 2.3 x 10 46 ergs s" 1 . We inter- 
polated the SED between the observed photometric points 
from 10 9 Hz to 10 18 Hz and derived a bolometric luminos- 
ity of — 2.4 x 10 42 ergs s" 1 , corresponding to a bolometric 
correction of -6.2 for the 8.0 /im IRAC band. The AGN 
therefore has a low Eddington ratio of — 10~ 4 . We also cal- 
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TABLE 1 

GALFIT DETERMINED PARAMETERS FOR SERSIC MODEL AND NUCLEAR POINT SOURCE 



Band 


Sersic Integrated Flux 


Sky 


Nuclear Flux 


Nuclear Luminosity 


(pm) 


(mjy) 


DC Offset 


(mjy) 


(10 41 ergss" 1 ) 


3.6 


2390 ±70 


-0.10 


4.67 ±0.22 


2.40 ±0.40 


4.5 


1420 ± 40 


0.57 


5.63 ±0.26 


2.32±0.38 


5.8 


1190 ±40 


-0.05 


6.78 ±0.27 


2.16±0.35 


8.0 


560 ±17 


0.11 


16.6±0.5 


3.85 ±0.62 


24 


400 ±12 


-0.20 


60.8±1.8 


4.69 ±0.76 



NOTE. — The Sersic model index (n= 6.07 ± 0. 10), effective radius (r eff = 146"2 ± 0''6), axis ratio (0.688 ± 0.004), and position angle (234.0° ±0.2° east of north) were all fitted 
only at 3.6 /j.m and were held constant for the fits made in the other bands. The central position of the Sersic profiles and the position of the nucleus were determined iteratively using 
the 3.6 fim and 8 jum images to be (03''22"'41 s .69, -37° 12'28"8) and (03 ,1 22'"4r ! .71. -37° 12'28"7) respectively. Note that the Sersic index determined matches well with the ones 
obtained by Cote et al. (2007). Luminosities in each photometric band given here are uL u in ergs s _1 , determined from the luminosity density L v . The uncertainty in the luminosities 
is the result of both the uncertainty in the flux measurement and in the distance to NGC 1316. 
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FIG. 3.— 3'0 X 3'0 images of NGC 1316: (a) Hubble ACS visible emission (535 nm), (b) Spitzer non-stellar 8.0 ^m emission, (c) VIA 4.89 GHz emission 
with a resolution of 1"38 X 1"02, (d) Chandra soft X-ray (0.5-2.0 keV) emission, (e) XMM-Newton soft X-ray (0.5-2.0 keV) emission, (f) Chandra hard X-ray 
(2.0-7.0 keV) emission. The Chandra images are smoothed with a 5" Gaussian and have pixel sizes of 1" and 0"5 for the soft and hard images respectively. The 
XMM-Newton image is the central region of Figuref4] demonstrating the north-south elongation of soft X-ray emission. With the exception of panel a, the darker 
regions have more emission. Panels a, c, and e have a logarithmic color scale, while panels b, d, and f have a linear color scale. The scale corresponds to 30" 
(3.3 kpc). These images show the variety of morphology present in this galaxy. Optical dust extinction coincides with infrared dust emission, but the distribution 
of the dust is distinctly different from that of the hot X-ray emitting gas. 



culated the Bondi accretion rate to be 1.6 x 1 0~ 4 M Pl yr" 1 
dBondi|[l952l) . based on the black h ole mass dNowak et al.l 
120081) . the 0.77 keV gas temperature dlsobe et alj|2006l) . and 
a central gas density of 0.4 cm" 3 derived from a /3 model 
dCavaliere & Fu sco-Femiano 1976) fit to the central 200" 
(21.7 kpc) region using Chandra data (ax = 320 pc, f3 = 0.49, 
no = 0.4 cm" 3 ) dJones et alj|20!ob . 

4.2. Inner Jet 



The 4.89 GHz radio emission on small scales (Fig- 
ure [3fc) has previous l y bee n extensively described by 
Geldzahler & Fomalontl(l l984l) . The northwest jet extends 0.'5 
(3.3 kpc) from the nucleus and does not decrease significantly 
in brightness over the initial 15" (1.6 kpc). In contrast, the 
weaker southeast counter] et decreases immediately in bright- 
ness away from the nucle us. Unlike the jet in M87 (IShi et al.1 
l2007llForman et alJ 2007) which is clearly detected in all four 
IRAC bands, the NGC 1316 radio jet is not detected in any 
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3CL0 24:00.0 3DJ0 3:23:DB.Q 3O.0 22:00.0 30JO 21:00.0 
- I I i , , , , I I . I ". 

FIG. 4. — 30' soft (0.5-2.0 keV) exposure-corrected XMM-Newton MOS image. Large X-ray cavities (shown by yellow circles) can be seen to the west (right) 
and to the southeast as lighter regions, centered at (3 h 22 m 16 s , -37° 14' 00") and (3''23 m 3\ -37° 15'40") respectively. The edges of the western cavity are marked 
by emission at (3 /! 22 m 28.5 s , -37° 17'42") and (3*22 m 24 s , -37°09'35"). The southeastern cavity has a faint edge to its north at (3*23"W, -37° 1 1'33"). These 
cavities contrast to the brighter emission to the north and south of the nucleus. The blue contours at 6, 11, 15, and 20 mjy/beam are from a 20 cm radio image 
and show the location of the radio lobes. The scale bar corresponds to 5' (32.6 kpc). 
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FIG. 5.— 3'0 X 3'0 images of NGC 1316: (a, b, c) Spitzer non-stellar 8.0 ^m emission with (a) contours of CO(2-l) (green) IHorellou et"aTll200Tb . (b) flat- 
fielded soft (0.5-2.0 keV) Chandra X-ray contours (blue) at 0.01, 0.02, 0.1, and 0.2 counts/ks/arcsec 2 , and (c) radio contours (red) at 0.04, 0.06, 0.10, 0.20, 0.50, 
1.0, and 6.0 mJy/beam (rms is 0.03 mJy/beam); and (d) Chandra soft X-ray (0.5-2.0 keV) emission with radio contours. The scale corresponds to 30" (3.3 kpc). 
These images show that: 1 ) infrared dust and molecular emission coincide, 2) the soft X-ray emission shows no indication of absorption due to the dust and cold 
gas, and 3) that the resolved radio jet does not coincide with dust or soft X-ray emission, but lies in an inner X-ray cavity. 
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TABLE 2 
Jet Infrared Flux Upper Limits 



Band 


Flux 


1 a Unc. 


Luminosity 


la Unc. 


(iim) 


(mJy) 


(mJy) 


(10 40 ergss" 1 ) 


(10 40 ergss" 1 ) 


3.6 


<1.05 


0.35 


<5.94 


1.98 


4.5 


<0.66 


0.22 


<2.98 


0.99 


5.8 


<0.58 


0.19 


<2.07 


0.69 


8.0 


<0.39 


0.13 


<0.97 


0.32 


24 


<0.81 


0.27 


<0.69 


0.23 



NOTE. — Fluxes, given in mJy, are 3cr upper limits and were obtained using an aper- 
ture denned by the radio emission northwest of the nucleus applied to each Spitzer im- 
age. The galaxy flux was removed using the Sersic model determined from fitting the 
Spitzer images with GALFIT. Luminosities in each photometric band given here are 
vL v in ergs s~ 3 , determined from the luminosity density L„ . 

IRAC band or at 24 /im. We used an aperture defined by the 
region of radio emission to the northwest of the nucleus to de- 
rive upper limits on the IR emission from the jet, which are 
given in Table [2] Extended aperture corrections were derived 
and appliedQ We also calculated the expected fluxes assum- 
ing a synchrotron model with a typical spectral index of 0.55 
from the radio flux of 29 mJy within the aperture. The ex- 
pected fluxes range from 0.38 mJy to 0.14 mJy, which are 
below the limits measurable with the present data. 

The soft X-ray image (Figure [3}l and [3^ and contours on 
Figure \5&) does not show emission from the jet. The re- 
solved northwest radio jet coincides with a region of low X- 
ray emission, most likely a small ~ 15" X-ray cavity created 
by the expansion of the r adio plasma, previously described by 
iKim & Fabbianol ( 12003b . As illustrated in FigureEJ;, the dust 
emission is faint at the position of the radio jet. The bend in 
the northwestern jet is located just south of the first IR knot, 
along the northwestern dust protrusion. 

4.3. Extended Non-Stellar Emission 

iTemi et aTl (120051) found that the morphology of the 
8.0 /im non-stellar emission was similar to that of the 
15 /im emission detected b y the Infrared Space Observatory 
(ISO) (Xil ouris et al.l2004l) . They concluded that, while much 
of the excess emission at 8.0 /im was likely due to PAH emis- 
sion at 7.7 /im, warm, small dust grains also contributed. The 
similarity of the features at 5.8 /im and 24 /im supports this 
interpretation. The extended non-stellar emission has signif- 
icant structure at all wavelengths as shown in Figure [3] and 
Figure [5] 

• The Hubble ACS image (Figure [3^) demonstrates that 
regions of visible dust extinction have similar morphol- 
ogy to the infrared dust emission described in §3.1 and 
shown at 8.0 um in Figure (See similar image in 
ITemi etalJd200l .1 

• The CO contours dHorellou et al.ll200lT) superimposed 
on the 8.0 /im non-stellar emission in Figure|5ji demon- 
strate that the northwestern and southeastern dust emis- 
sion regions coincide with molecular hydrogen traced 
by the CO emission, suggesting a common origin for 
the dust and cold gas. 

• The strongest X-ray emission outside the nucleus 
(Figure |5j) with soft X-ray contours overlaid on the 

6 IR AC: Extended Source Calibration 

http://ssc.spitzer.caltech.edu/irac/iracinstrumenthandbook/33 



8.0 /im non-stellar emission) extends northeast of the 
nucleus along the major axis of NGC 1316 into a region 
absent of infrared dust emission. We tested whether the 
dusty features seen in the infrared and the coincident 
cold gas result in soft X-ray absorption, but found no 
indication thereof, although the data is not sufficient to 
make a conclusive statement. 

In the galaxy core, the soft Chandra image and the XMM- 
Newton image show a roughly north-south elongation approx- 
imately 1 .'25 (8.2 kpc) in each direction (Figures [3}i, |3£), 
which does not follow the distribution of the stars. Instead, 
this emission is roughly perpendicular to the major axis of 
NGC 1316 and may be from hot gas that was moved by 
the outburst. The larger XMM-Newton field of view (Figure 
|4| shows further filamentary emission north of the nucleus 
and a pair of X-ray cavities. These cavities, likely created 
by the expansion of radio plasma, are marked with yellow 
circles of 230" (25 kpc) radii. The western cavity is cen- 
tered at (3 A 22 m 16 s , -37° WOO") and the southeastern cavity 
is centered at (3*23 m 3 s , -37°15'40"). Each cavity lies 320" 
(34.8 kpc in the plane of the sky) from the nucleus. There are 
three regions of enhanced emission along the edges of these 
cavities, approximately located at (3 A 22 m 28.5 s , -37°17'42"), 
(3 ft 22 m 24*, -37°09'35"), and (3 /, 23 m 00 I , -37°11'33"), which 
are likely due to increased gas density as the hot ISM is com- 
pressed by the expanding cavities. No radio emission is de- 
tected in the se X-ray cavities, a situati on previously seen i n 
Abell 4059 dHeinz et all l2002h. M87 (iForman et all l200l . 
and the Perseus cluster (Fabian et al. 2006). While the centers 
of the radio lobes line up with the AGN, there are indications 
that this system may be experiencing some sloshing of the hot 
gas. Specifically, the X-r ay cavities are cen tered \ '.5 and 3 '2 
south of the nucleus and Eke rs et al.l (119831) found low-level 
radio emission between the lobes ~ 7' south of the nucleus. 

To quantitatively measure the significance of the cavities 
seen in the XMM-Newton image, we plot in Figure [7] the az- 
imuthal surface brightness of the soft (0.5-2.0 keV) exposure- 
corrected background-subtracted XMM-Newton image (Fig- 
ure HJ taken in an annulus extending from 180" (19.6 kpc) 
to 375" (40.8 kpc) from the nucleus after the image was 
smoothed with a 28. "2 Gaussian. The azimuthal profile shows 
that the regions between 80-140°, (i.e. southeast of the nu- 
cleus), and between 220-280°, (i.e. to the west), have signif- 
icantly lower surface brightness than the rest of the annulus. 
These regions coincide with the cavities identified (yellow cir- 
cles) in Figure H] The bright regions north of the nucleus and 
along the cavity edges in Figure |4] coincide with the signifi- 
cantly brighter regions in the azimuthal plot (Figure [7]). The 
three regions of enhanced X-ray emission along the cavity 
edges are noted in Figure [7] 

We tested whether the variations in azimuthal surface 
brightness could be the result of abundance or gas density 
variations. The maximum surface brightness change would 
require a factor of 2.3 difference in elemental abundance (i.e. 
the lower surface brightness region would have an elemental 
abundance 40% that of the brighter regions). While such an 
abundance gradient would be relatively long lived against dif- 
fusion, even if it proceed s as fast as predicted for hea vy ions 
in a fully ionized plasma (Sarazin 1988, Spitzer 1956), such a 
distribution of metals mimicking cavity structures seems par- 
ticularly contrived. An alternative explanation for the surface 
brightness variations is for the isobaric gas to have a density in 
the regions of lower surface brightness 0.66 times that of the 
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FIG. 6. — Spectral energy distribution of the nuclear emission in the Spitzer and Chandra bandpasses, through the F814W Hubble filter (upper limit), at 4.9 GHz 
and 15.0 GHz from Geldzahler & Fomalont 1 1984), and at 1730A from Fabbiano et al. ( 1994). The right vertical axis applies to this AGN, and the left axis to 
the comparison AGN SEDs that are normalized at 1 fim and come from Figure 7 of Ho (2008). The purple squares are the LLAGN with \og(Ltoi /LecIcI < —3.0). 
NGC 1316's AGN has log(L io ; /L Edd = -3.9). Note the large radio to optical ratio and the lack of the UV-optical big blue bump in the NGC 1316 AGN. 
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FIG. 7. — Azimuthal profile of background-subtracted counts per square 
arcsecond for aradial sector with inner and outer radii of 180" (19.6 kpc) and 
375" (40.8 kpc) from the nucleus in the soft (0.5-2.0 keV) XMM-Newton im- 
age in Figure [4] Angles are given counterclockwise from north. Note that 
there is less flux to the southeast (120°) and to the west (250°), which cor- 
responds to the locations of the centers of the cavities apparent in Figure [4] 



gas to the north and south of the nucleus. Such a difference 
in density requires either that the lower surface brightness gas 
is at least 1.5 times hotter than the surrounding gas or that 
the regions of lower surface brightness would be filled with 
a relativistic plasma. No indication of emission from a hot 
plasma is seen in the harder X-ray band (2.0-7.0 keV). There- 
fore, we expect the regions of lower surface brightness to be 



cavities filled with a currently undetected relativistic plasma, 
a morphology seen in oth er galaxies and ga laxy clusters and 
known as ghost cavities (Hei nz et all 120021) . Emission from 
ghost cavities are generally detected in low-frequency radio 
data (e.g. Giacintucci et al. 2009). 

5. MULTIWAVELENGTH VIEW OF THE 
CONSEQUENCES OF THE MERGER EVENT 

NGC 1316 exhibits signs of a recent merger, including nu- 
clear activity and a disturbed morphology seen in the opti- 
cal and i nfrared dus t distri bution as well as the tidal tails first 
noted bv lSchweizerl {1980). Each wavelength provides differ- 
ent insights into the merger event and the resulting structure of 
NGC 1316. Below, we discuss the distribution of the infrared- 
emitting dust and estimate the mass of the galaxy that collided 
withNGC 1316from the measured dustmass. We also use the 
morphology of the large scale radio and X-ray emission (Fig- 
ure @|l to constrain the recent outburst history of the central 
SMBH. 

5.1. Dust Distribution 

To measure the dust emission in the Spitzer bands, we per- 
formed aperture photometry on the infrared images after sub- 
tracting a Sersic model of the stellar emission. Elliptical aper- 
tures (Figure [8) were chosen to include dust features seen 
at 8.0 /im. For the southeastern region, we used an ellipse 
with a major axis of 56."8, a minor axis of 46."6, and a posi- 
tion angle of 284° (east of north) centered at (3 /! 22 m 42 s .75, - 
37° 13'02"1). For the northwestern region, we used an ellipse 
with a major axis of 86", a minor axis of 54", and a posi- 
tion angle of 245° centered at (3''22"W.47, -37°1 1'53."0). 
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FIG. 8. — 8.0 fim non-stellar image showing the position and size of the 
apertures used to determine the dust photometry. The southeastern aperture 
is an ellipse with a major axis of 56" 8, a minor axis of 46"6, and a position 
angle of 284° (east of north) centered at (3*22 m 42 s .75, -37° 13' 02'.' 1). The 
northwestern aperture is an ellipse with a major axis of 86", a minor axis of 
54", and a position angle of 245°centered at (3*22 m 40 1 .47, -37° 1 1'53"0). 
The scale bar corresponds to 30" (3.3 kpc). 



TABLE 3 
Dust Photometry 



Band 


SE Knot 


NW Knot+Arc 


Total 


Total 




03*22"*42 s .75 


03*22'"40 1 .47 


of 


with 




-37°13'02?'l 


-37°11'53"0 


knots 


nucleus 1 


(/im) 


(mjy) 


(mJy) 


(mjy) 


(mJy) 


5.8 


4.16±0.43 


6.47±0.70 


10.6±0.8 


17.4±0.9 


8.0 


18.8±0.6 


13.8±0.5 


32.6±0.8 


49.2±0.9 



NOTE. — Fluxes (cols 2-3) were obtained using elliptical apertures (Figure [8j for 
the southeast and northwest region on each Spitzer image once a Sersic model had been 
subtracted. 

' Nuclear fluxes are given in TablefTI 



The counts in each aperture were background subtracted and 
converted to fluxes, after an extended aperture correction was 
applied. Columns 2-3 of Table [3] list the flux in each aper- 
ture. The next two columns give the total fluxes and their un- 
certainties with and without the nuclear point source (nuclear 
fluxes given in Table [TJ. The uncertainties have two compo- 
nents: the uncertainty in the Sersic model and th e uncertainty 
in the photometric accuracy of Spitzer images (Reach et alj 
120051 E ngelbracht et 10 120071) . 

IDraine et all (12007b modeled the integrated IRAC, MIPS, 
and IRAS fluxes of the SINGS galaxy sample with a two- 
component dust model and determined dust masses. They es- 
timate their models are accurate to 10%. To test whether the 
dust mass determined for NGC 1316 by IDraine et al.l (120071) 
is located within the regions visible in Figure[8] we compared 
our photometry from the two regions of dust emission and the 
nucleus t o their dust model p redictions. We convolved the 
predicted IDraine et al.l d2007l) dust flux SED for NGC 1316 
(their Figure 14) through the appropriate response functions to 
calculate the expected fluxes in the Spitzer band passes. Our 



5.8 /im and 8.0 /im total fluxes for the dust regions and nu- 
cle us of 17.4±0 . 9 mJy and 49.2±0.9 mJy are consistent with 
the IDraine elaT] d2007l) model values of 13.0 ± 1.3 mJy and 
45.4 ±4.5 mJy, which are the emission for the entire galaxy. 
The agre ement leads us to co nclude that the dust mass esti- 
mated by IDraine et al.l (120071) is contained in the non-stellar 
IR emission regions described in §3.1 and in the nuclear re- 
gion. 

5.2. Dust Mass and Merger Progenitor Mass 

In the following, we show that the dust observed in NGC 
1316 is not native to the galaxy and use the dust mass to esti- 
mate the mass of the merger galaxy. The clumpy morphology 
of the dust is significantly different from the smooth ellipti- 
cal distribution of the stars, so the NGC 1316 s t ars co uld not 
have expelled the dust. In addition, iTang et aT] (|2009) found 
for nearby ellipticals that most of the non-stellar 8.0 /im emis- 
sion is confined to the nuclear region. This also demonstrates 
that th e morphology of th e dust emission in NGC 1316 is un- 
usual. Te mi et alJ ([2009) found a correlation between the K- 
bandand24 /im luminosities of elliptical galaxies. NGC 1316 
has a particularly large 24 /im luminosity for its K-band lumi- 
nos ity, about an order of magnitude greate r than predi c ted by 
the iTemi et all (l2009h correlation. While iTemi et all (120091) 
did not find a correlation between K-band luminosity and ei- 
ther 70 /im or 160 /im luminosity, NGC 1316's integrated lu- 
minosities at these wavelengths of 1 .4 x 10 43 ergs s -1 and 
1.5 x 10 43 ergs s" 1 dDale et al.ll2007l) are also more than an 
order of magnitu de greater than found for the galaxies in the 
ITemi et ail (120091) sample. The large infrared luminosities of 
NGC 1316 demonstrate an external origin for the dusty emis- 
sion. 

In the following, we estimate the mass of dust in NGC 
1316 as well as the dust mass e xpected to be in an el l iptica l 
galaxy the size of NGC 1316. iMufioz-Mateos et alj (2009) 
provided a formula (A8) for calculating the dust mass of a 
galaxy from its 24 /im, 70 /im, and 160 /im fluxes and its 
distance. We calculated that NGC 1316 has a total dust mass 
of 2.4 ±0 9 x 10 7 M , using the integrated MIPS fluxes from 
iDaleet all (120071) . While emission from the dust is clearly ob- 
served, the large uncertainty on the dust mass results from the 
uncertainties in the distance to NGC 1316 (22.7 ± 1.8 Mpc) 
and in the integrated MIPS fluxes (0.43 ± 0.02 Jy at 24 fjxn, 
5.44 ±0.40 Jy at 70 /im, and 12 .61 ± 1.78 Jy at 160 /mi). 
We revised the IDraine etai] J2007) dust mass for NGC 1316, 
which was found on the basis of SED fitting, for our assumed 
distance of 22.7 Mpc to be 3.2 x 10 7 M Q , which is consis- 
tent with our dust ma ss. We used the fl uxes for the sample 
of elliptica l galaxies inlTemi et al.l (120091) . along with their B- 
V colors (Ide Vaucouleurs et alJ 19911) an d the color-dependent 
mass-to-light ratios of lBell et aD Q2003), to calculate the stel- 
lar and dust masses of the sample. We found that ellipti- 
cal galaxies typically have dust-to-stellar mass ratios between 
0.7-5.3 x 10" 7 . Using these ratios, we estimate that NGC 
1316 wi th its stellar mass of 5.3 x 10 11 M Q ( based on B- 
V=0. 87 dde Vaucouleurs et aljfl99ll). K = 5.5 87 dJarrett et alj 
120031) . and the relations of iBeiTeT al. (2003)) had an intrin- 
sic dust mass of 0.4-3 x 10 5 M Q , < 1% of the measured dust 
mass. We conclude that nearly all of the dust currently present 
in NGC 1316 was contributed by a merger galaxy. 

We can constrain the galaxy type and the stellar and gas 
mass of the merger galaxy from its estimated dust mass of 
2.4 ±0.9 x 10 7 M . The merger galaxy had to be a late 
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type galaxy as its stellar mass, were it a typical elliptical, 
would have been roughly 200 times the present stellar mass 
of NGC 1316. We calculated the stellar masses of the spi- 
ral galaxies in the SINGS sample u sing NED colors and the 
color-dependent stell ar M/L ratios oflBell et alJd2003l) . Using 
the dust masses from Draine et al . (2007), we found dust-to- 
stellar mass ratios between 0.4-3.4 x 10~ 3 for Sa-Sm galaxies 
with Sc galaxies having the largest ratios. From these ratios 
and the merger galaxy dust mass, we estimate that the merger 
galaxy had a stellar mass in the range of 1 - 6 x 10 10 M Q , 
approximately 10% of NGC 1 316's current s tellar mass. As- 
suming typical galaxy colors (Trimble 2000), we calculated 
the co rresponding Lr t o be 0.7-2 x 10 10 L Q using the ra- 
tios of Bell et al] (120031) . We estimate the corresp onding cold 
gas masses based on gas-mass-to-light ratios of Betto niet al.l 
d2003|)tobe2-4x 10 9 . iKennicutt et all d2003l) found an 
upper limit on the mass of neutral hydrogen in NGC 1316 of 
5.5 x 10 8 M , and lHorellouet all (1200 ll) estimate that NGC 
1316 has 7.4 x 10 8 M Q of molecular hydrogen gas, resulting 
in a total cold gas mass of less than 1.3 x 10 9 M Q . Since the 
merger galaxy's estimated cold gas mass of 2-4 x 10 9 M Q is 
larger than NGC 1316's present cold gas mass, some cold gas 
may have been ionized due to mixing with hot gas or used in 
star formation in the merger process. 

5.3. Outburst Ages 

The morphologies of the large scale GHz radio emission 
and the large X-ray cavities in NGC 1316 suggest two possi- 
ble interpretations of the recent outburst history of the NGC 
1316 AGN. Based on the relative location of the X-ray cav- 
ities and the radio lobes, we can conclude that either the 
1.4 GHz radio features do not fully define the extent of the 
radio lobes created in conjunction with the X-ray cavities in 
the course of a single outburst or there were at least two out- 
bursts, one resulting in the radio lobes and a more recent one 
creating the X-ray cavities seen in the XMM-Newton image. 
In the two outburst scenario, we expect the X-ray cavities, 
which lie at a smaller radius, to result from the more recent 
outburst, as they would otherwise have been disrupted by the 
expanding radio lobes. 

Cavity ages can be estimated by assuming that the bub- 
bles that c reate them rise buoya ntly in the gaseous atmo- 
sphere (e.g. lChurazov et ai]|2001l) . Assuming an approximate 
buoyancy velocity of ~0.6 c s , where c s is the sound speed 
dChurazov et al.ll2001l) . we estimate a buoyancy speed of 270 
km s" 1 in the 0.77 keV medium (llsobe et al.l2006j). The X-ray 
cavities are located 320" from the nucleus and lWadd (11961b 
measured the separation of the radio lobes to be 33'. If we 
assume the lobe expansion is in the plane of the sky, for NGC 
1316's distance of 22.7 Mpc, the X-ray cavities and the ra- 
dio lobes are at 35 kpc and 108 kpc from the nucleus respec- 
tively. These distances correspond to buoyancy rise times of 
0.1 Gyr for the X-ray cavities and 0.4 Gyr for the radio lobes. 
The age of 0.4 Gyr appears large for radio lobes still emit- 
ting at 1 .4 GHz, b ut matches the es timate of the synchrotron 
age calculated by lEkers et al.l dl983l) . Further, if the lobe is 
continuously or intermittently connected to the nuclear power 
supply, fresh injection of electrons or re-energization of ex- 
isting electron populations could, in principle, supply suffi- 
cient high energy electrons so that the lobe is visible at GHz 
frequencies even afte r several hundred million years. Finally, 
llvomoto et al. ( 1 998 ) estimated that the nucleus of NGC 1316 
was active ~0. 1 Gyr ago, which agrees with our estimate of 



th e age of the X-ray cavities. 

Mackie & Fabbianoj (1 1998b estimated a single merger with 
a low-mass gas-rich galaxy ^0.5 Gyr ago could be respon- 
sible for the optical tidal tail morphology. Such a merger 
could have provided the material to power the AGN outbursts. 
A gas-rich merger galaxy would contribute blue young stars, 
which require about a Gyr to become red and dead. The B- 
V color of NGC 1316 of 0.87 dde Vaucouleurs et alJfl99ll) is 
slightly blue r than the typ ical B-V color of 0.91 for ellipti- 
cal galaxies dTrimble II2000T) , which suggests that NGC 1316 
may still contain a small population of early type stars. NGC 
1316's color therefore supports a merger within the last Gyr. 
We also can set a lower limit on the merger age by estimat- 
ing the free -fall time of the northwe stern dust component. We 
revised the Arnabol di et al.l dJ998) total mass within 45" of 
the nucleus to be 6.6 x 10 10 Mq and used this along with 
the excess velocit y of 70 km s" 1 for th e northwestern clump 
of molecular gas dHorellou et a l. 2001) to estimate a free-fall 
time for the northwestern dust and molecular gas feature. The 
estimated free-fall time of 22 Myr is a lower limit on the age 
of the merger, since the clump likely also has a tangential ve- 
locity component and was likely deposited at greater radii by 
the merger. 

5.4. Outburst Energies 

iChurazov et al.l d2002l) described the energy deposition re- 
quired to inflate a bubble adiabatically, and thereby create a 
cavity, as the enthalpy of that bubble, which for relativistic 
gas is 4PV . To estimate the energy of the outburst responsible 
for the X-ray cavities, we used the more clearly defined west- 
ern X-ray cavity, whose shape we approximate as a sphere of 
230" (25 kpc) centered 320" (34.8 kpc) from the nucleus. To 
meas ure the pressure, we assume an isothermal gas at 0.77 
keV (jlsobe et al. 2006) and solar abundance and derived the 
density from the surface brightness. We model the density 
as a /3 model (Cavaliere & Fusco-Femiano 1976), whose pa- 
rameters we derive by fitting the surface brightness profile of 
the exposure-corrected XMM-Newton image in a region not 
containing the cavities. This method provides a lower limit 
on the total outburst energy, since it estimates the kinetic en- 
ergy released in the outburst. Assuming solar abundance, we 
estimate the kinetic outburst energy is 10 58 ergs, for equal- 
sized bubbles created in the plane of the sky on each side of 
the nucleus. If the abundance is half solar, then the gas den- 
sity and outburst energy both increase by ^40 %. Based on 
the energy needed to create the X-ray cavities and adopting 
a mass-energy conversion efficiency of e = 10%, we estimate 
the mass of material that would have been accreted onto the 
SMBH to be: 

AM BH =- -^ = 5xlO 4 M (1) 

e c z 

where E is the total energy output. 

Deep X-ray observations are not available for a similar 
analysis in the regions defined by the radio lobes, and we ex- 
pect complications due to inverse Co mpton X-ray emissio n 
coincident with the potential cavities (Feigelson et alJI 1995b . 
However, if we extrapolate the gas density model to the ra- 
dio lobes assuming they lie in the plane of the sky and use 
the 20 cm observation to determine the location and size of 
the lobes, we can estimate the energy required to evacuate 
cavities the present size of the lobes and thereby estimate the 
energy of the outburst required to create them. We approxi- 
mated the lobes as 24' (78.3 kpc) diameter spheres centered 
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at 14.'3 (93.3 kpc) (west) and 15.'6 (101.4 kpc) (east) from 
the nucleus. We extrapolated the gas density model derived 
from the XMM-Newton surface brightness to the radii of the 
radio lobes and combined the derived pressure there with the 
expected cavity volumes created through adiabatic bubble ex- 
pansion to estimate the required energy of ~ 5 x 10 58 ergs. 
The mass accreted onto the SMBH to produce this energy 
would be ~2 x 10 5 M Q . 

5.5. Comparison of the Cen A and Fornax A Jets 

Centaurus A (NGC 5128) and NGC 1316 (Fornax A) are 
both nearby elliptical galaxies, which have recently under- 
gone a merger event that has produced strong nuclear ac- 
tivity. They each host a 10 s M© black hole that are low 
luminosity AGNs and ha ve dust lanes rough l y perpendic- 
ular t o their radio lobes ( Ma rconi et al.l 120061 iNowak et al.l 
2008). However, they differ significantly with regards to the 
observatio nal characteristics o f their jets. At a distance of 
3.7 Mpc dFerrarese et al.l 120071) . Cen A has a 1.5 kpc radio 
and X-ray emitting j et that extends from the nucleus to the 
northeast radio lobe dFeigelson et al.lll981l iKraft et al.|[2002i 
lHardcastle et al.ll2003l) . In contrast. NGC 1316's ~3 kpc ra- 
dio jet does not extend to the large radio lobes and coincides 
not with an X-ray jet, but with a soft X-ray cavity. This phe- 
nomenologica l contrast suggests that, while the jet of Cen A 
is dissipative dNulsen et al]|2010l) . NGC 1316's jet is not or 
that the NGC 1316 jet has shut off on larger scales. 

6. CONCLUSIONS 

We detected considerably more dust emission for NGC 
1316 than expected in an early type galaxy with its K-band 
luminosity and observed evidence of recent AGN outbursts 
in the form of X-ray cavities and radio lobes. We presented 
Spitzer images of the infrared dust emission, including the 
first image of dust emission at 5.8 /xm. We determined that 
the dust has a clumpy morphology, mostly confined to two 
regions, one 28 "8 (3.1 kpc) southeast of the nucleus and a 
43 "9 (4.8 kpc) protrusion ending in an arc northwest of the 
nucleus. Molecular emission is detected from these regions. 
The resolved radio jet is not detected by Spitzer, and it does 
not coincide with regions of dust emission. Since the dust 
must be almost entirely external in origin, we use the dust 
mass to constrain the type and mass of the merger galaxy. We 
calculated a present dust mass of 2.4 ±0.9 x 10 7 M© based 
on the integrated MIPS fluxes, which agrees with the dust 



mass o f 3.2 x 10 7 M Q predicted by the model of Drai ne et"afl 
d2007l) . We estimate the merger galaxy was a late type galaxy 
with a stellar mass of 1-6 x 10 10 M© and a gas mass of 
2-4 x 10 9 M©, some of which was likely ionized or used 
to form stars in the merger event. 

We constrained the age and energy of the merger and out- 
burst events based on the X-ray and radio emission. The 
XMM-Newton image shows a pair of X-ray cavities to the 
west and southeast of the nucleus, likely created by the ex- 
pansion of radio plasma, which are closer to the nucleus than 
the 1.4 GHz radio lobes. The relative locations of these cavi- 
ties and the radio lobes suggests that either the 1 .4 GHz radio 
emission does not show the full extent of the radio emission 
from the outburst or that there have been at least two AGN 
outbursts. We calculate buoyant rise times for the X-ray cav- 
ities of 0. 1 Gyr and for the radio lobes of 0.4 Gyr, assuming 
expansion in the plane of the sky at 0.6 c s , which agr ees with 
the synchrotron age estimated by Eker s et al.l (119831) . Since 
the age of the radio lobes is close to the 0.5 Gyr age esti- 
mated by Mackie & Fab bianol d!998l) for the merger, the out- 
burst was likely triggered by the accretion of material onto 
the SMBH from this merger. Finally, we constrained the ki- 
netic energy of the outbursts based on the energy required to 
create the XMM-Newton cavities and the radio lobes bubbles. 
We estimate the outburst that created the X-ray cavities had a 
kinetic energy of 10 58 ergs and that the creation of the radio 
lobes required ~ 5 times more power. 
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